Revisiting Bremsstrahlung emission associated with Light Dark Matter annihilations 
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We compute the single bremsstrahlung emission associated with the pair annihilation of spin-0 particles 
into electrons and positrons, via the r-channel exchange of a heavy fermion. We compare our result with the 
work of Beacom et al. Ql- Unlike what is stated in the literature, we show that the Bremsstrahlung cross 
section is not necessarily given by the tree-level annihilation cross section (for a generalized kinematics) times 
a factor related to the emission of a soft photon. Such a factorization appears only in the soft photon limit or 
in the approximation where the masses of the particles in the initial and final states are negligible with respect 
to the mass of the internal particle. However, in the latter case, we do not recover the same factor as for 
e+e^ fj'^fj^j. Numerically the difference, in the hard photon limit, is as large as a factor 3.6. However the 
effect on the upper limit of the dark matter mass is not significant. Using gamma ray observations, we obtain an 
upper limit on the dark matter mass of 30 or 100 MeV (depending on the region that is considered for the 511 
keV analysis) while Beacom et al. found 20 MeV. 



1. INTRODUCTION AND SET UP 

The Milky Way radiation has been mapped extensively at dif- 
ferent wavelengths for many decades. Photon emission at low 
and ultra-high energy is still currently studied. All the latest 
observations appear extremely exciting and shed new light on 
the physics of the galactic centre. Among them, there is the 
observation of a 5 1 1 keV emission line by SPI, a spectrometer 
on board of the INTEGRAL satellite launched in 2002. The 
511 ke V line has been detected by many experiments since 
the seventies (see Ref. 10, IH 01), but the origin of the line 
could not be established unambiguously. With the new results 
from SPI |5, 6, 7], one could determine that it is due to (para- 
)positronium formation, thus confirming the presence of anti 
matter in our galaxy. 

A very striking feature of this line is that it is mostly emit- 
ted in the bulge of the Milky Way (the distribution is well 
approximated by a sphere; the full width at half maximum is 
about 10° of diameter). No significant emission was detected 
in the disk. This characteristic tends to exclude astrophysi- 
cal sources such as hypernovaeAVolf-Rayet stars, pulsars, and 
cosmic ray interactions. On the other hand, this observation 
could point towards an old galactic population jsl]. 
A blind source analysis showed that more than 8 point sources 
could explain the emission 0|. A diffuse source provides a 
very good fit too. None of these sources has been found as 
yet. This makes the origin of the galactic positrons even more 
mysterious. 

Low Mass X-ray Binaries (LMXB) could perhaps be an ex- 
planation if the positrons emitted in the disk escape into the 
bulge. Their characteristics are poorly known though (e.g. 
flux, jet content). No 511 keV emission has been seen so 
far from the brightest of these objects, but more observa- 
tions will be dedicated to this kind of candidates. Type la 
supernovae (SNla) may also explain the line emission if the 
positron escape fraction and the supernovae explosion rate are 
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large enough to maintain a steady flux. At present various ob- 
servations contradict with each other and it is quite difficult to 
know whether these two conditions are indeed satisfied in our 
galaxy |39]. Dedicated observations of SNla remnants are 
also scheduled with INTEGRAL and should help answering 
these questions l^lstl- 

Another possible source of low energy positrons could be dark 
matter annihilations 1 10]. This scenario requires a dark matter 
mass between a few MeV and a few hundred MeV, two kinds 
of interactions, and a cuspy dark matter halo profile (likely 
to be as cuspy as a Navarro-Frenk- White model 1 11]). This 
dark matter profile coiTesponds to theoretical predictions from 
numerical simulations, but it has not been confirmed as yet 
observationally ll o.ri3ll . 

Light dark matter (LDM) particles were initially proposed in 
the context of the so-called cold dark matter crisis. These 
particles were meant to experience a new damping effect, in- 
cluding one that is half-way between the collisional damp- 
ing and the free-streaming. They were meant to inherit from 
the damping of relativistic species (X = y,v) with which they 
weakly interact either directly or indirectly fl4l ITsl fl^ flTl 
[1^ 1. In the indirect case, the damping of the species X is com- 
municated to dark matter through the dark matter coupling to 
other particles in contact with X. For example, before recom- 
bination, dark matter indirectly experiences the Silk damping, 
until it thermally decouples from electrons. This new effect 
(called induced damping) is simply the generalization of the 
Silk damping effect to dark matter and has been confirmed nu- 
merically ] 19]. In the direct case, LDM may suffer from the 
mixed damping effect ] 14, 15, 17], which corresponds to the 
situation where, due to a weak but large enough neutrino-DM 
interaction rate, dark matter inherits from the neutrino free- 
streaming. This occurs when LDM stays coupled to neutrinos 
after their thermal decoupling. For some specific values of the 
DM-v elastic scattering cross section, these interactions affect 
the DM fluid but do not change the neutrino fluid properties. 
This effect was stressed in Refs. ] 14, 15] for neutralinos when 
lliey are strongly mass degenerated with sneutrinos and for 
MeV particles having weak interactions n^fT^fTsil . 
Light dark matter particles therefore have this interesting 
property: on the one hand, they look like cold dark matter 
candidates because of their mass and weak interactions. On 
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the other hand, their linear matter power spectrum is cut-off 
at a cosmological scale (due to these new damping effects), 
which is a well-known characteristic of warm dark matter can- 
didates. 

To avoid an overproduction of gamma-rays and satisfy simul- 
taneously the relic density criterion, the LDM total annihila- 
tion cross section at early times should be larger than that in 
the primordial universe or in virialized objects 1 20] . This leads 
to a scenario in which the LDM interactions, and hence the 
damping, are suppressed fTs^ (the damping mass in the linear 
matter power spectrum was estimated to be about 100 Mq). 
Following up Ref. fioll . there have been many other dark mat- 
ter models proposed to explain the 511 keV emission, e.g. 
Refs. |21, 22, 23, 24, 25, 26]. Here we focus on the annihilat- 
ing case. In this scenario, there is a light extra gauge boson Z' 
(needed for the relic de nsity ) and heavy charged particles (for 
the 511 keV emission) 12011 . This new set of particles is quite 
similar to that expected in = 2 supersymmetry. 
The model is described in detail in Re f. It27iil. The constraints 
on the couplings are given in Ref. iHiilsl. We disregard 
the possibility that dark matter is fermionic I40ll . following up 
Ref. [M- 

Heavy dark matter particles, i.e. 100 MeV > mdm > 10 
MeV, are in danger of yielding a too large contribution to the 
anomalous magnetic moment of the electron ]29, 30]. Also 
they would make the observed fraction of para and ortho- 
positronium very hard to explain and would eventually 
yield a too large e+e^ inflight annihilation contribution ]31]. 
Any other way of reducing the mass range which is not too 
model dependent and which could be related to the 511 ke V 
flux is very interesting. In that context, the Beacom et al. pro- 
posal 1 1 ], to constrain the dark matter mass from the compari- 
son between observations and the gamma-ray flux originating 
from the Bremsstrahlung, looks promising even if it gives a 
smaller contribution than the inflight annihilations. 
This suggestion is complementary to a previous estimate 
made before INTEGRAL'S publication ]20]. In that work, it 
was assumed that each electron and positron produced by dark 
matter annihilations spontaneously generates one photon after 
their emission. The authors assumed that the total annihila- 
tion cross section was either velocity-dependent or velocity- 
independent. After comparison with the gamma-ray observa- 
tions (and requiring the correct relic density), they concluded 
that only a velocity-dependent cross section was possible for 
Wdm < 100 MeV while for greater masses either a velocity- 
dependent or independent cross section could be allowed. 
However, the model that fits INTEGRAL/SPI data requires 
two kinds of interactions simultaneously. One is velocity- 
dependent and is needed for fitting the relic density. The other 
one is velocity-independent and is quite mandatory for fitting 
the 511 keV emission morphology. Such a scheme is not com- 
patible with Ref. ]20] in which it was assumed that only one 
process (either velocity dependent or independent) was con- 
tributing to the total annihilation cross section. One therefore 
has to estimate the gamma ray production again and derive 
the corresponding upper bound on the dark matter mass, as- 
suming that LDM is at the origin of all the 511 keV emission 
mapped by SPI. 



By doing so the authors in Ref. fl*] concluded that mdm could 
not exceed 20 MeV. Their estimate of a^^ drn'^p+^-y b^sed 
on an old observation made by Berends et al. ]32] according 
to which, in QED or SU{N) gauge theories, the cross sections 
associated with the single emission of a hard bremsstrahlung 
photon becomes remarkably simple in the ultrarelativistic 
limit. Using specific examples, Berends et al. showed that 
a factorization into two terms can be observed. One factor is 
the squared amplitude associated with the 2-to-2 process (i.e. 
without the additional emission of the photon), and evaluated 
for a generalized kinematics. The second factor is related to 
the infrared factor, which describes the soft photon limit. Such 
a property shows up for example in e+e^ /u^/u^j. The ex- 
pressions for the initial and final-state radiation cross sections 
associated with this process can be found for example in the 
review of Berends and Bohm ] 33| . 

Beacom et al. thus argued that the final-state radiation cross 
section of dark matter particles annihilating while they are al- 
most at rest can be written as the product of the two-to-two 
process times an universal factor. In a next step they assumed 
that this factor can be obtained from the known results for the 
final-state radiation associated with the e+e^ — > /j+/j^y pro- 
cess in the ultrarelativistic limit, where the electron mass is 
neglected. 

However, in the case of dark matter, the ultrarelativistic limit 
is not appropriate. Furthermore for chiral dark matter cou- 
plings it is also obvious that the electron mass cannot be ne- 
glected. Also e+e^ ("^/J ^"d dm dm* e+e^ are intrinsi- 
cally different, e+e^ — > is due to a photon exchange in 
the i-channel; it involves only vectorial couplings. In contrast, 
dm dm* e+e^ proceeds through a f-channel exchange of a 
heavy fermion and can involve chirality flipping couplings. 
Thus writing the Bremsstrahlung cross section as the product 
of two factors may not be correct. 

Here we demonstrate that the approach followed in Ref. (l|] 
fails to correctly describe the Bremsstrahlung process in the 
hard photon limit. 

II. ANNIHILATION AND BREMSSTRAHLUNG CROSS 
SECTION 




FIG. 1 : Feynman diagrams contribution to 5 + 5* — > e+ + e + Y- 

In what follows we study the process 

S{pdm) + S* (/5dm* {p,- ) + e+ ) + y(/Jy) , ( 1 ) 
where S denotes a spin-0 dark matter particle. 
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Since dark matter is neutral, the emission of a photon neces- 
sarily comes from the final state particles and the particle that 
is exchanged during the annihilation process (F). The corre- 
sponding Feynman diagrams are shown in Fig. ^ It is conve- 
nient to use the scaled energies as independent variables: 



and 



2k -p^- 



2k ■ 



2k- 



Py 



(2) 



where k = p^m + Pdm* = Pe- + Pe+ + Py and s = (pdm + 
Pdm*)^- Due to energy momentum conservation, x,x, and Xy 
are not independent of each other: 

2=x + x + Xy. (3) 

In terms of x,x and Xy, the scalar products between final state 
momenta read: 



Pe+ -Py 




(4) 


Pe- -Py = 




(5) 


^e--Pe+ = 


-{1-Xy-2Z), 


(6) 



where z — ntg^/s and is the electron mass. In the case of 
dark matter it is legitimate to assume that it is almost at rest 
when it annihilates. At Earth position, the dark matter veloc- 
ity dispersion is about three orders of magnitude smaller than 
the speed of light (v ~ lO^^c). It is even smaller below 1 kpc 
(say, inside the galactic centre). Thus any scalar product that 
involves the four- momenta of both initial and final state parti- 
cles can be expressed in a very simple way, without involving 
any undefined angles. In particular we obtain: 



(7) 
(8) 
(9) 



A. The annihilation cross section 

For later use we reproduce here the annihilation cross section 
associated with the process 
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(10) 



at tree-level 12711 . Neglecting the P-wave annihilation cross 
section because of the suppression of the dark matter velocity 
in the 511 ke V region, we find that the corresponding cross 
section Oo is given by: 
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with 



(11) 



(12) 



4m, 



dm 



(13) 



where m^^ denotes the mass of the annihilating scalar and nip 
the mass of the heavy charged fermion that is exchanged in the 
f -channel. From inspection of Eq. il 1> . one can easily see that 
in general it is not appropriate to neglect the electron mass. 
If the ratio between the chiral couplings c/ , is of the order 
of nie/mf then the term that is formally suppressed through 
the electron mass has the same importance as the term in nif. 
Such a scenario arises naturally in 2-Higgs doublet models 
with flavour changing neutral couplings. I.e. the lagrangian 

in ffif 

L = -gi—eRFLS-gr—FReLS + h.c. (14) 

V V 

leads immediately to ^ = ^ if we identify ci=gi'-^ and ~ 
gr'-^- Such a case does not fit SPI data unless one considers 
extremely large couplings of a few units. Since the particle F 
is charged, the mass mp must be larger than 100 GeV, given 
LEP limits 1 34] restricting the range of the couplings that can 
be considered. 

In the limit where and <C nif, the cross section be- 
comes: 



Oo = 



:^^-^^{chCr,me,mi^,mF) + o( , (15) 
3271 p,- niF \mF^ J 



with 



K{ci,Cr,me,mdm,mF) 



2ciCr+{ci^+Cr^) 



me \ ' 

niF ) 



2 ime'^-mAm^ 
8Cr Q 5 . (16) 

niF" 



Keeping only the leading term in the l/niF expansion, we ob- 
tain for Oo 



oVPi^^M 



(17) 



To compute the annihilation rate and therefore the flux of 
positrons and/or photons that are emitted during the dark mat- 
ter annihilations, one must consider Oov^ (where is the dark 
matter relative velocity) instead of Oq. v, is twice the dark 
matter velocity in the centre-of-mass frame (written here as 
P,). We therefore recover the results of Ref. Ii27il . 



B. Bremsstrahlung cross section 

1. General formulae 

In terms of the squared matrix element My, the single differ- 
ential bremsstrahlung cross section reads 



dOy 

dXy 



1 



2567t3p; 



dx |Myp, 



(18) 
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where the bounds x_ ,x+ are given by 



(19) 



The squared matrix element is obtained from the evaluation of 
the Feynman diagrams shown in Fig. Since there are many 
different mass scales, the complete result - keeping the full 
mass dependence - is rather lengthy. We find: 



2 '-m/ 



(iciCrMF + (c/^ + Cr^ 
-8ci'^Cr^(^{2xy- l)mdn 

+ SctCr{ci^ + c/)memF^'S 

— 4cic,.{ci^ + c,-^)mi,mF(^m/ +mim^{2xy— l)jc 
^ (20) 



where e denotes the electric charge and where the explicit re- 
sults for the functions J?,!B,C,©,£,5\t are given in the ap- 
pendix. We note that in addition to the symmetry in c;,Cr, 
there is also a symmetry inx^ x, when using x instead of Xy. 
Furthermore we have checked that we reproduce the correct 
factorization when the emitted photon becomes soft: 



\My 



iPe- ■ Pi) {Pe+ ■ Pi) (Pe- ' Py)^ {Pe+ ' Pyf 



\Mo\\ 
(21) 



For arbitrary hard photons no such factorization can in general 
be observed. 

There are many extra terms in the squared amplitude with re- 
spect to that of the annihilation. Some are particularly relevant 
when nif /2 < mdm ^ nif (e.g. XyA, 'D, iAt) or in the case 
of very large couplings, i.e. {cj + c^)me ^ 2ciCrmf. 
Hence from Eq. i20i one readily sees that - unless one makes 
simplifying assumptions - the cross section for the emission 
of an additional photon cannot be written as the tree-level 
cross section for a generalized kinematics times a factor re- 
lated to the soft photon emission. 

The integration of the squared matrix element \My\^ over x to 
obtain the photon spectrum is straightforward. In the follow- 
ing we only present analytic results in the approximation that 
the dark matter mass is small with respect to the mass of the 
particle F. 



2. Results in the limit mp 2> (mdm, nig) 

Expanding Eq.|5U|in terms of 1 /nif and keeping the dominant 
terms up to we obtain: 

K{ci,Cr,s,me,mF) 



My 



1 



1 



{l-x-xyf (l-x)2 
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(22) 



Integrating over x we obtain the photon spectrum: 
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(23) 



Keeping the electron mass only in the logarithm so as to regu- 
larize the mass singularity, the above result can be simplified 
further: 



da 



'J 



dXy K X' 



ao^lUl + ^)lnfA)-2^^ (24) 



where we introduced 

■s' = (Pe- +Pe+)^=S (1 -Xy). 



(25) 



For e+ + e ^ /j+ + /J + y, keeping only the final-state ra- 
diation, the formula is given by 1 33]: 



do 



dXy 



= a, 



n'2 



(26) 



Comparing Eq. J24> and Eq. i26\ . it is easy to see that the 
bremsstrahlung cross section associated with the annihilation 
of two scalars cannot be obtained from the final-state radiation 
in muon pair production. However, in the soft photon limit 
(i' s) the two results agree after the identification 
nie, as it must be. We are thus led to the conclusion that the 
bremsstrahlung correction to the annihilation of two scalars 
cannot be evaluated as it was done in Ref. 



III. DISCUSSION AND UPPER LIMIT ON THE DARK 
MATTER MASS 

In the upper plot in Figure |2l we show the Bremsstrahlung 
cross section as given in Ref. 1 1 ] (dashed lines) versus our ex- 
pression (dotted lines) for 2, 20, 30 and 1 00 MeV. All our plots 
are normalized to Oo ^ . To show more explicitly the differ- 
ence between these two expressions, we show in the lower 
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plot the ratio: 



R = 



^^-''WEq.|24| 



for 1, 10 and 100 MeV. 
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in Ref. f)S- 

obtain: 



For maximal photon energy (xy = 1 — 4z), we 



(l-4z)2 



(1 + I6z2) (ln(4)-l)' 



(27) 



In the limit z 0, Eg. 1271 tends to 3.6. Although the differ- 
ence between Ref. Ql] and the results presented in this let- 
ter becomes larger for large mdm> the phenomenological rel- 
evance might be larger for smaller m^i^. For large m^i^, the 
ratio is peaked at the phase-space boundary in contrast with 
small where a larger phase space region is affected. 
The upper limit on the dark matter mass comes from the 
hard photon region. Fig. |3l shows the expected flux of 
Bremsstrahlung photons for dark matter particles with a mass 
of 20, 30 and 100 MeV and for < 5°, |Z| < 5° (dashed lines) 
and \b\ < 5°, |/| < 30° (soUd Hnes) with b the latitude and / 
the longitude. The expected flux in Ref. 1 1] is represented by 
the dotted Unes for mdm = 30 MeV. 




10 
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FIG. 3: Bremsstrahlung flux for 20, 30 and 100 MeV, in the limit 
(OTdm ^ riif ), for b < 5° and / < 5° (dashed lines) and b < 5° and 
/ < 30° (solid lines). The data are represented by diamonds and the 
result from Ref. QJ is plotted for comparison (dotted lines). 



FIG. 2: Upper plot: our Bremsstrahlung cross section (dotted lines) 
versus the Bremsstrahlung cross section used in Ref. (Jl (dashed 
lines) for 2, 20, 30 and 100 MeV. Lower plot: ratio of the two ex- 
pressions for 1, 10 and 100 MeV. 

The lower plot in Figure |2] shows that the result obtained in 
this paper is up to a factor 3.6 larger than the result obtained 



Despite the sizeable difference in the hard photon limit the up- 
per limit on the dark matter mass is not changed significantly 
due to the limited precision of the data. We obtain the same 
result as in |30], i.e. < 30 or 100 MeV, depending on 
whether one uses SPI data from regions defined by \b\ < 5° 
and |/| < 5° (for which the 511 keV line intensity is about 
0.018 ph/cm^/s) or |/| < 30° (for which the intensity is 0.005 
ph/cm^/s) respectively. 
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Many papers, e.g. Refs. fss'.'s^,'??] use the expression that is 
given in 1 1]. It might be worth checking whether their results 
remain unchanged using the exact formula presented in this 
work. 



IV. CONCLUSIONS 

In light of Ref. we computed the Bremsstrahlung emission 
associated with the annihilation of spin-0 particle dark matter 
particles and compared it with the result for e+e^ /j+/j^y. 
We find that, in general, one cannot write the Bremsstrahlung 
cross section as the tree-level annihilation cross section times 
a factor related to soft photon emission. In the limit where the 
mass of the new charged fermion F (which is exchanged in 
the f-channel) is much heavier than the dark matter mass mdm 
and the electron mass Wg, we observe a factorization similar 
to that assumed by Beacom et al. However, the factor we 
find from the explicit calculation is different from that given 
in Ref. yj. We thus conclude that the bremsstrahlung cross 
section associated with the annihilation of spin-0 particles into 
an electron-positron pair cannot be obtained from final-state 
radiation in e+ + ^ /j+ + /j^ . 

The difference between Ref. and our result is sizeable. 
However, the analysis presented in Ref. f3^ remains valid. 
For b <5° and / < 5°, we find that the upper limit on the dark 
matter mass is about 30 MeV. For b < 5° and / < 30°, it is 
rather 100 MeV. 
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Appendix 

The squared amplitude for the Bremsstrahlung process is 
given by: 



|Mf|2 



Di 



_^ ^ Re(M,-M;^) ^ ^ Re(MfM;_) ^ ^Re(MfM;^) 



DfD, 



(28) 



where M; ~ ee*^ u^^Mi Vg and 

^e- = 'Y^,(Jy + ^e-+me){clP| + CrPR) 

X (Je-+^y-^dm + mF){ciPR + CrPl), 

M,+ = {c,Pi + CrPR) (j>,- - /dm + mp) 

X {c,Pr + CrPl ) ( -/y - + me)yf, , 

•Mf = (c/P/ + CrPR) {J>e- -/dm + mp) 

X liA {^Am* - + fflF ) {ciPr + CrPl), 



(29) 

(30) 
(31) 



with 
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S = Mg Xy — 2mdm m (1 — x-).)(2xx — 2 + 2ji:y+x-). ) 
+ mdm^me^ [ - 8 + 24xy - 2(}Xy + 9xy 



xi(-8 + 16Xy- llXy^)] 



^dm'[-4/(l 



-%x\2-l,Xy 



C = ntg X. 



(-20 + 36xy - 1 8x/ + ) 

Xy(-4+ 10Xy-9Xy2 + 3Xy^) 
x(8-12Xy+8Xy^-3Xy'*)], 
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Y 



y ~V Ay 



— 2OTdm^'We'^(l ^Xy) [ — 2 — 2x2 + 2x(2 — Xy) + 2x- 
+ 2mdm''(l -x){l -Xy)[-2 + 4x^ + 2xy + Xy2-Xy^ 
+ 4x2(-3 + 2xy)+x(10- 12xy+3xy2)] 

4 2 r 2 3 

+ mdm fng [ — 8 + 24xy— 27xy + lOxy 

- 4x^(2 - 4xy+ 3xy2) -4x(-4 + lOxy- 10xy2 + 3Xy^)] , 

— nig^ll ~ 2Xy)Xy^ 

+ 2m(]m^(l — 2xy)2(l — X — Xy)(2 — 2Xy + Xy2) 

+ 2OTdm^'We'^ [2 — 8xy + 8xy2 — 2xy^ + x^ (2 — 6xy + 3xy2 ) 

+ x(-4 + 14Xy- 12Xy2 + 3Xy^)] 

+ mdm'^'We^ [ ^ 16 - 8x'*(l - Xy) + 64xy- 97xy2 + 68xy^ 

- 18xy'^+ 16x^(2 -3xy + Xy2)+2x2(-28 + 64xy-49xy2 



15x7 



2x(24 - 76xy+ 90xy'^ - 5 Ixy^ 



llXy^)], 



2mdm'*[(l — 2x)miai^ + nig^ — niF^]^ {I — x)^ 
[(1 — 2x)OTdm^ + mg^ — mF^]^{l — x)2, 
—2ci^Cr^mg^[mg^ + mdm^{l — 2x)] 
[m/ + mdm^ {l-2x)]^ + (c,'* + c/) 1^ 



%, = OTe%2 + 2mdm'*(l --«)(! --=c)[2 + Xy(-6 + 5xy)] 
+ ni^m^fng^ [4 — 4xx(l — Xy) +Xy(— 8 +Xy(3 + 2xy))] , 



1^ = 



— OTf, '"Xy2 — 2OTdm**'W£'^ (1 ^X)(l — 2Xy) 



X [ — 6 + 8xx(l — Xy) + (10 — 3Xy)Xy] 



- 16mdm"^(l -xf{\ -xf{\ -xy)[(l -x)2 + (1 -xf 

-4 + 4xx(l -Xy) -Xy(-8+Xy+4Xy2)] 
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mdm mg 
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12-2xx(6+xy(-14+ llXy)) 



7 



xy ( 40 +Xy( - 47 + 2x^(7 + 2xy)) 

A-\6x^{-l+Xy)- 32x^ (2 - xy) [l-Xy) 
+ JCy ( 8 + Xy(-5 1 + 8 (8 - 3xy)xy) 



I 6 4 



4x^( 21 +jCY(-37+XY(14 + 3xy)) 



4X(2 - Xy) 5 +Xy(-5 +Xy(-6 + 7Xy)) 



(32) 
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